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ABSTRACT 

We examine the effects of streaming cosmic rays upstream of a strong, parallel collisionless shock. 
We include explicitly the inertia of the cosmic rays in our analysis, which was neglected in previous 
work. For parameters relevant to the acceleration of cosmic rays at a supernova blast wave, we find no 
MHD fluid instability that would lead to the amplification of the magnetic field above that given by 
the compression at the shock. We show how to recover, from our own analysis, the cosmic-ray-driven 
MHD fluid instability found by previous authors. We conclude that including the inertia and dynamics 
of the cosmic rays keeps the system stable. More over, the cosmic ray current leads to an additional 
Hall-like term in the magnetic evolution equation. The implications of this paper for acceleration of 
galactic cosmic rays at supernova remnants are briefly discussed. 
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1. INTRODUCTION 

Blast waves propagating out from supernova remnants 
(hereinafter SNR) are thought to be the source of most 
galactic cosmic rays (hereinafter GCR) up to an energy 
of ~ 10 15 eV (IKoyama"eFaT1 Il995t JAllen et all Il997t 
Tanimori et al.lll998l IAharonianlll999t iBerezhko fc Volkl 
2003t IVink and Laming] 120031 ). Because the observed 
GCR spectrum follows a power law with constant index 
up to this energy, diffusive shock acceleration (hereinafter 
PSA) is the most likely primary acceleration mechanism 
(lAxford et all U977I: [Belli [l978t iBlandford and Ostrikerl 
119781: Uokipiil U982L I1987Q . For acceleration in a SNR, 
DSA at a quasi-parallel shock is only able to account 
for energies up to ~ 10 14 eV using the observed in- 
terstellar magnetic field a nd maximum scattering rate 
(jLagage &; Cesarsky||1983l) . In order to reach the knee 
at 10 15 eV the acceleration occurring at such a shock 
must be more rapid. In the case of DSA, this means the 
particles must be bound to the shock more effectively by 
the magnetic field. In the case of a parallel shock, more 
rapid acceleration can be accomplished by increasing the 
magni t ude of the magnetic field. 

iBelll (|2004f ) has suggested a magnetohydrodynamic 
(hereinafter MHD) fluid instability which amplifies 
a transverse perturbation upstream of the shock to 
SB /Bo ~ 100, where Bo is the magnitude of the original 
field. This cosmic ray current-driven instability (here- 
inafter CRCD) is triggered by GCRs streaming along the 
uniform background magnetic field. This GCR current 
produces a separate return current in the background 
plasma which excites transverse turbulence in the up- 
stream region. 

The sole purpose of the GCRs is to generate the re- 
turn current which maintains charge neutrality. The 
GCR inertia is not considered as we will explain be- 
low. Once established, the excited transverse compo- 
nent scatters high energy GCRs, limiting their motion 
away from the shock increasing their acceleration rate. 
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He finds that the amplified field is adequate to acceler- 
ate GCRs up to, and possibly beyond, 10 15 eV. Obser- 
vations of supernova remnants (Bcrezhko & Volk I2003T) 
strongly support the existence of an amplified mag- 
netic field of the required magnitude. Others have also 
worked in the kinetic regime to examine Bell's stream- 
ing instability. While some work supports the large 
increase m SB (jReville et all 120061: iBlasi and Amatol 
2008), some find the amplification is only moderate 
SB / Bp < 10 (iRiauelme and Spitkovsky| |2009h. Recently, 
iRiquelme and Spitkovskv ( 2010) have proposed an addi- 
tional instability caused by the perpendicular streaming 
of GCRs relative to the background field. This perpen- 
dicular current-driven instability (hereinafter PDCI) is 
actually triggered by the pre-amplified field generated 
by the CRCD. Additional work has also been carried 
out concerning k inetic inst abilities related to self excited 
waves m shocks (|Ledll983D . 

An additional possibility to account for the observed 
acceleration to 10 15 eV is that much of the acceleration 
takes place where the shock is quasi-perpendicular, in 
which case no upstream magn etic-field amplification is 
necessary (| Jokipiil fl98i LT987I) . In that case, the ob- 
served amplification of the magnetic field could be the 
result of preexisting upstream turbulence warping the 
shock front, leading to turbulent motions of the plasma 
behind the shock that ampli fy the magnetic field through 
a turbulent dynamo pr ocess (jGiacalone and Jokipiill2007l : 
iBeresnvak et al.ll2009l) . 

In the present paper we will take a new approach to 
the original MHD problem of streaming GCRs at a quasi- 
parallel shock and compare the results with those found 
previously. We find that if the inertia of the streaming 
cosmic rays is included and their dynamics considered, 
there is no MHD fluid instability. 



2. LINEAR ANALYSIS 



2.1. 



The CRCD MHD Fluid Instability 

It is well known that energetic particles accelerated at 
a quasi-parallel shock can stream ahead of the shock and 
excite waves in the upstream medium through kinetic 
processes. Recently it has been proposed that these same 
particles can excite MHD fluid instabilities in the back- 
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ground fluid as well. In this case, the shock-accelerated 
GCRs stream along the assumed uniform background 
magnetic field with a bulk speed v s , the speed of the 
shock responsible for initially accelerating the GCRs. As 
mentioned above this results in the amplification of mag- 
netic perturbations transverse to the background mag- 
netic field. 

Beginning with the momentum e quati on for the b ack- 
ground fluid. ILucek and Belli (l2000h and|]pi3 (l200l find 
that the GCRs drive a MHD fluid instability in the back- 
ground fluid. Note that in this section and in direct refer- 
ences to lBelll ((20043, results will be presented in SI units 
as in the original work. 



du 



-VP+j x B 



(1) 



Here p, P and B are the mass density, pressure and mag- 
netic field of the background fluid. The background fluid 
moves with velocity u and the current in the background 
fluid is j. 

The streaming GCRs are now introduced by way of 
Ampere's law3, VxB = /io(j + jcr — n cr eu) where the 
RHS of this expression is the total current in the system. 
Here j cr is the current due to the streaming GCRs and 
the third term on the RHS preserves the quasi-neutrality 
of the system by providing an additional population of 
background electrons which cancel the positive charge 
of the GCRs. The excess electrons are made to move 
with the same velocity as the background fluid, u. Sub- 
stituting for the background current j using the above 
definition of the total current neglects the GCR inertia 
while including the forces acting on them and leads to the 
modified momentum equation for the background fluid. 



du 

} ~dt 



VP 
1 

Po 



( V x B) x B - j cr x B + n cr eu x B (2) 



The second line of equation ^ is equal to the magnetic 
Lorentz force on the background plasma after solving 
for j in Ampere's law above. The term — j cr x B will 
ultimately drive the instability as discussed below. At 
this point, we stress that the inertial term on the LHS of 
equation (|2|) represents only the inertia of the background 
fluid with mass density p and velocity u. A correspond- 
ing inertial term for the GCRs (with variables p cr and 
u cr ) is not present in this equation and only the electric 
current due to the GCRs has been retained. This GCR 
current produces the return current in the background 
plasma which drives the instability. In this way, the in- 
ertia of the GCRs has been neglected in Bell's treatment. 
All subsequent discussions of inertia relate to this term. 

In order to time-evolve the magnetic field in Bell's 
treatment, the assumption is made that the GCRs have 
no effect on the magnetic induction equation and the 
magnetic field is frozen into the background plasma as in 
ideal MHD. 



^Vx(uxB), 



(3) 



3 Note the use [iq (due to use of SI units) in the above definition 
of the current. 



Equations ^ and © are then linearized using the ex- 
pressions B = B|| + Bj^, u = ui, E = Ej_ and 
jcr = jii +j_L- All values labeled || are time-independent 
zcroth order terms oriented in the z direction (taken to 
be the direction of the uniform magnetic field) and those 
labeled _L are first order perturbations in the x — y plane 
which vary in time and space. The gradient of the pres- 
sure is neglected at this point because it does not affect 
the linear study of the problem at hand. Equations ^ 
and ^ are then combined to eliminate leaving a dif- 
ferential equation for the perturb ed m a gnetic field and 
perturbed current (equation (3) in[Bell (2004])). 



d 2 Bj 

dt 2 



,d 2 Bj 

1 dz 2 



"JII 



PV S 

<9B_l 

dz 



B 



<9B, 



Of 



p dz 



x B|| = 



(4) 



The speeds v a and v s are the background Alfven and 
shock speeds, respectively. The perturbed current is 
eliminated using the relationship jj^ = (aj i|/B|i)B_l 
where a is a uni tless quantity defined below (equation 
(4) inlBel l|200"l ). 




v±Pa 



kvn 



(5) 



dp\\ — 

u df 
k dp\ 



Here e, p, v, and / are the GCR charge, GCR momen- 
tum, GCR velocity and the GCR di stribution f unction, 
respectively. This is the term used by Bell (2004) to cou- 
ple the GCRs to the magnetic field though this effect 
will be neglected in the study of the CRCD instability 
(discussed in section 13.1 p . As before, _L and 1 1 denote 
quantities perpendicular and parallel to the equilibrium 
magnetic field. 

Assuming an oscillatory solution for Bj_ in equation 
(@]), ~ exp(ikz — iujt), B ell finds the following dispersion 
relation (equation (7) in Bell (2004)): 



u 2 -v 2 a k 2 ±(v 2 s — (1^0x0 = 



(6) 



Here, £ is a unitless parameter quantifying the strength 
with which the GCRs drive the instability. 



c = 



j\\ B \\ r ai 
pv 2 s 



(7) 



This variable follows from the j cr x B term in ©. The 
subscript I on r g i, the GCR gyro radius, and a refer 
to those values calculated for GCRs with momentum pi, 
the lower bound on the pre-existing CR spectrum. The 
quantity o\ is related to but differ^ from a above. 

The consequences of equation (|6]) will be discussed in 
section 13.11 

2.2. Composite Fluid 

In our present treatment of this problem, all of the in- 
dividual particle populations are treated as fluids which 

4 See lEeH [12001) Sections 3 and 4 for details 
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are then combined into a single fluid. From this point 
forward, we refer to this single fluid as a composite fluid. 
The inertia of the GCRs (here assumed to be singly 
charged with rest mass m p ) is explicitly included by con- 
sidering the dynamical equation of the GCRs. This is in 
contrast to the work of previous authors working in the 
MHD limit who have ignored the inertia of the GCRs 
and simply included the compensating return current in 
the momentum equation of the background fluid (as in 
equation ([2]) above) . We begin with four separate particle 
populations which are the background ions (n b i), back- 
ground electrons {n\, e = n^), GCRs (n cr ) and the excess 
electrons required for quasi- neutrality (n ex = n cr ). The 
separation of the excess electrons from those in the back- 
ground is not necessary but done here to emphasize the 
quasi-neutrality of the composite fluid. 

2.2.1. Composite Momentum Equation 

The composite momentum equation is found by adding 
together the momentum equations of the four popula- 
tions. The index a implies a sum over the four popula- 
tions. 



4 d 4 (\ 

^ —n a m a u a = ^2 n ae a ( -u Q x B + E 



a = l 



(8) 



For simplicity, we have taken the total plasma pressure 
to be isotropic. Carrying out the summation gives the 
following composite momentum equation. 



du h 



n bl m p - 



dt 



du be du. 
ntrim e —^- + (7 cr ) n cr m p 



dt 



dt 



d\i P 



- (n bi u bi ~ n bi u be 
c 



TL cr \l n 



c) xB 



dt 



(9) 



The derivative used above is the total or convective 
derivative, d/dt a = d/dt + u a ■ V where u a is the ve- 
locity of the fluid element being followed. The velocity 
subscripts hi, foe, cr, and ex stand for background ions, 
background electrons, cosmics rays and excess electrons 
respectively. The RHS of this equation is the total cur- 
rent, j, crossed with the magnetic field. The factor of 
{"far) in the above expressions arises due to the use of 
the relativistic momentum for the GCRs and represents 
an average Lorentz factor for the population of GCRs. 
While the bulk speed of the GCRs, u cr , is not relativis- 
tic, individual GCRs are relativistic in the shock frame 
while being accelerated by DSA. The GCR fluid is es- 
sentially treated as a fluid of heavy ions enhancing their 
contribution to the momentum of the composite fluid. 
Our work is not particularly sensitive to this quantity so 
we will simply treat (7 cr ) as a constant and vary its value 
over a wide range of values. A more in depth description 
of (7cr) is outlined in Appendix [Bl 

Summing the le ft han d side as outlined in 
iBovd and Sanderson! (|1969ft (Appendix [A) equation 
J9J simplifies to: 



du 1 ._ 
^ = 4^ (VXB)XB 



(10) 



Here u is the center of mass velocity and p is the total 
mass density of the composite fluid. Note that the term 
on the LHS of equation (|TU|) includes the inertia of the 



GCRs. Upon neglecting the contributions of the elec- 
trons to the fluid momentum (rn e /rrii << 1), u and p 
can be expressed as below. 



u 



nuxi b i 



cr j Tlcr^-cr 



n bi 



<7cr) n c 



P =(n bi + (7 cr ) n cr )m p 



-Tl r n7Tlj j 



(11) 



(12) 



We have substituted for the total current j using Am- 
pere's law in the non-relativistic MHD approximation. 



—V x B 

47T 



(13) 



Including the dynamical equation of the GCRs from 
the beginning we avoid the more complicated momen- 
tum equation of Bell in equation ((2} where the forces on 
the GCRs are included without their associated inertia. 
Instead we find the usual momentum equation from ideal 
MHD. 

Due to the approximation of isotropic pressure the fire- 
hose instability will not appear in our analysis. It should 
also be not ed that the kinetic resonan t streaming GCR 
instability (jKulsrud and Pearcel 1 19691 ) is also absent in 
our r esults due to our use of the M HD fluid equations. 
[Belli ([200l and lReville eTaTl ([200(1 both note, however, 
that the growth of the non-resonant instability dominates 
that of the resonant instability. 

2.2.2. Composite Induction Equation 

Rather than assuming that the GCRs have no ef- 
fect on the MHD induction equation, we explicitly de- 
termine the MHD electric field of this composite fluid 
(|Bovd and Sanderson! [2001 . To do so, we begin with 
the MHD momentum equation for the electrons. 

m e ^-(n bi u be + n cr u ex ) (14) 
dt 

= --{nbi\i be + n cr u ex ) x B e(n bi + n cr )E 
c 

Now, in the limit that m e « m, we neglect the the 
inertial term on the left hand side of equation (fl"5|) . 

« - {n bi u be + n cr u ex ) x B - e(n bi + n cr )E (15) 
c 

Using the above equation we can solve for the electric 
field, E, to find a generalized Ohm's Law. 



E 



1 

n e c ' 



) x B 



(16) 



This shows the magnetic field is actually frozen into the 
total electron fluid as in the case of Hall MHD. We define 
the variable n e = n b i + n cr . We then eliminate the elec- 
tron velocities using the definition of the total current, 
j = e(n bl u bl + n cr u cr - n bl u be - n cr u ex ). 



ln bi 1 . 
E = um x B H j x B 

c n e n e ec 



1 n r 



-u cr x B (17) 



By doing so we gain what we call an effective Hall term 
(the final two terms in equations (|17|) V The middle term 
is due to the transformation from the electron frame of 
reference to that of the ions and GCRs (center of mass) 
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and is present in Hall MHD. The final term on the RHS of 
(fl~7|) is due to the presence of the GCRs. Physically, the 
Hall terms represents an electric field that is generated 
when particles of opposite charge separated due to the 
Lorentz force in a purely magnetic field; the Hall effect. 
It can be thought of as a correction to the induction 
equation because the magnetic field is actually frozen 
into the electron fluid rather than the composite fluid 
moving at the center of mass velocity. The correction 
accounts for the ions and GCRs having much larger gyro- 
radii than the electrons. Lastly we express E in terms of 
the center of mass velocity u. 



E 



1 

n e c 
1 



[n m u + (1 - (Tor))" 

B 

(V x B) x B 



47m e e 

Inserting this expression for E into Faraday's law: 



(18) 



<9B 

~dt 



cV x E 



(19) 



we arrive at the following induction equation for the evo- 
lution of the magnetic field. 



dB ( n m 

9F =Vx ITT uxB 



- V x 
+ V x 



47rn e e 

(1 - (lcr))n 



(V x B) x B 

u rr x B 



(20) 



This differs from Bell's induction equation in © where 
the magnetic field is assumed to be frozen into the back- 
ground fluid. There, the GCRs were assumed to have 
no effect on the form of the induction equation. Equa- 
tion (|20p does not assume the magnetic field lines are 
frozen into the composite fluid and includes the effects 
of the GCRs. The second and third terms on the RHS of 
equation (|20[) prevent the composite fluid from being well 
magnetized. The only assumption necessary to arrive at 
this result is that m e /m p << 1. Consideration of the 
relativistic mass of the GCRs leads to the second term 
in (f20|) because their contribution to the momentum of 
the fluid differs from their contribution to the current by 
a factor of (jcr)- 

2.2.3. GCR Momentum Equation 

In order to find the perturbation to u cr (Su cr ) in our 
subsequent linear analysis of equations (fTU|) and ([2D]) we 
include the GCR momentum equation, (l2~Tj) below. This 
is necessary to provide a third equation for our system of 
three unknowns, SB, Su and Su cr . This method differs 
from Bell's calculation of o\ from (|5|) and was chosen as 
a more straight forward way to calculate the perturba- 
tions to the GCR velocity. Because Bell takes a% — >• 
to find his instability condition, the particular method 
for perturbing u cr is not critical to our comparison with 
Bell. What we consider to be important is retaining all 



effects due to the GCRs throughout our calculations. 



du cr ( 1 

(7cr) n cr m p —^- = n cr e I -u cr x B + E 



(21) 



The electric field can be eliminated by substituting E 
from equation (IT51) . 



/ v d\l cr 71m e . . 
(7cr) — TT = ( u cr - U) X B 



dt n e nipC 



4irn e m p 



(V x B) x B (22) 



2.2.4. Perturbation Analysts 

We now perturb equations ([TO]). and with 

transverse fluctuations in B, u and u cr . 



B = SB(it ± iy)exp(ikz — iu)t) + BqZ 



(23) 



u = Su(x ± iy)exp(ikz — iut) + (7 cr ) — —u cr z (24) 



<5ti cr (x ± iy)exp(ikz — iuji) + u cr z 



(25) 



Equations (f2"3"|) - ([2"5l) assume the initial magnetic field to 
be uniform in the z direction and the streaming GCRs 
(also along z) to be the only part of the composite fluid 
in motion. Both Bq and treated as constants. 

Because of the coupling of perturbations in the x and 
y directions we have assumed the perturbations to be 
circularly-polarized. The ± indicates the handedness of 
the perturbation with the + being right-handed and the 
— being left-handed. The wave number k is defined to be 
real while w can be complex. Both can take on positive 
or negative values. 

Inserting the assumed forms of B, u and u cr from 
above into equations (fTU)) , (I2U1) . and (|2"2")l we find the 
following system of equations in our three unknowns. 



nu kv 2 a 



SB 



U - (icr) ku c 



Su = 



(26) 



n cr Tt\ji k v 

oj ku cr =p 

71 e 



II r 



0J U 



6B + 


kB 




71 e 



Su 



n r 



(1 - { lcr ))^LkB 

Tie 



n b j kvl n bi LOgUcr 

Tie -B n e B 



5u rr = 



(27) 



(7c,.) UJ - (j cr ) ku c 



SBt 



i 

± OJg 

Tie 



Su 

Su cr = 



(28) 



In the above equations (|26|) follow from (|T0|. (1271) from 
([20)1 and ([28]) from ([22]). The variables v a and lo q are the 
Alfven speed in the background plasma and gyro period 
of the background ions respectively. Our assumption of 
transverse perturbations allows for the density parame- 
ters (n m and n e ) in equation (|20[) to remain constant in 
space and move outside the calculated curls. 

This system of equations can be expressed as a matrix 
equation, A ■ X = 0. 
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A solution for this system exists if the determinant of 
the coefficient matrix, A, vanishes. The Greek variables 
in A can be read from equations (|26|) - (l28|) . Calculating 
this determinant gives a dispersion relation that is cubic 
in uj. While the roots of a cubic equation have analytic 
solutions, the complexity of the polynomial coefficients 
results in complicated roots better suited to numerical 
methods. Rather than calculating the roots explicitly 
we will simply determine whether the roots will be purely 
real or complex over a wide range of the relevant param- 
eters k, jcr, u cr and n cr . To do so we calculated the 
variables Q and R. 



Q = 



36 



R 



2a 3 - 9ab + 27c 
54 



(29) 
(30) 



The variables a, b and c are the real coefficients of the 
monic cubic polynomial uj 3 +auj 2 +buj+c = 0. If the quan- 
tity Q 3 — R 2 < the roots of the cubic will be complex 
and h ave a non-zero imaginary component (j Press et al.l 
1992). In this case, the system would be unstable to 
perturbations. This condition will be evaluated over the 
wide range of parameters listed below. 



3 x l(T b G 

3 



Bo = 

n bl = 1 cm 
1 < (Tor) < 10 6 
lCP 15 rjf,i < n cr < lCF 1 ?!^ 

20^ < u cr < c/10 
-10 4 fc cr < k < 10 4 k cr 

Th e quantity k cr = 2ir /r g \ where r g \ — 1.1 x 10 18 cm 
from [Belli pOOl . 

In order to evaluate the condition for complex roots 
the values of {"f C r), n cr , u cr and k were varied over the 
ranges listed above by different step sizes. Both {^ cr ) 
and n cr were increased by powers of 10 and the value of 
u cr was increased by a factor of 20 with each successive 
calculation. Lastly, the value of k was incremented by 
10" 7 /e 

cr ■ 

In order to check the accuracy of our calculation of the 
condition Q 3 — R 2 < 0, (j C r) was set equal to 1 while 
all other parameters were varied according to the pro- 
cess outlined above. For this special case the roots of the 
dispersion relation are known to be real under all circum- 
stances so Q 3 — R 2 should never evaluate to a value < 0. 
Upon running the calculation we find several cases that 
return Q 3 — R 2 ~ — 10~ 31 , leading to false positives for 
complex roots. Even if these complex roots were valid, 
the resulting growth rates are on the order of 10 -12 s _1 
which are far too slow for the problem at hand. It would 
take 10 4 years to achieve <5B/B ~ 1. Using this in- 
formation we adjust the condition for complex roots to 
be Q 3 — R 2 < — 10~ 30 to account for numerical errors. 
Under this corrected condition we find all roots of the 
dispersion relation to be real. This means that there are 
no MHD fluid instabilities in our composite fluid formu- 
lation under reasonable conditions. 

3. COMPARISON WITH PREVIOUS WORK 



3.1. The Original Cosmic-Ray-Driven MHD Instability 

Bell (2004) considers three separate regimes in which 
to examine the dispersion relation of equation (|6|). He 
states that two of the regimes are of little interest as 
the waves are only weakly driven by the GCRs. These 
regimes are kr g i < 1 and kr g i > (^v 2 /v 2 . The intermedi- 
ate regime, 1 < kr g \ < (v 2 /v 2 , results in the GCR driven 
instability of interest here. In this case, his dispersion re- 
la tion s i mplifi es to uj 2 — v 2 k 2 ± (v 2 k/r g i (equation (15) 
of iBelll ((200 4)') under the approximation that o\ — > 0. 
The choice of + or — is related to polarization and de- 
termines whether or not the mode will be oscillatory or 
growing. A condition for the onset of the instability is 
readily found upon solving for lj in Bell's equation (15) 
above. We have substituted for £ in the equation below. 



uj = v a k. 



kB 



II 



(31) 



In order for the magnetic field to be unstable, the expres- 
sion under the radical must be negative to give a complex 
uj. This requires ct rninirnurci value for j cr . 



Jcr ^* 



Mo 



(32) 



Working under the assumption that the condition in 
is satisfied, Bell finds the maximum growth rate and cor 
responding wave number. 



In 



Cv 2 



2v a r„i 



Cv 2 



(33) 



(34) 



Recall that C represents the driving strength of the in- 
stability and follows from the j cr xB term in @ which 
only appears in ((2|) because the GCR current is intro- 
duced into the momentum equation of the background 
fluid without considering the GCR inertia. 

3.2. Simplifications of Composite Fluid to Recover 
Bell's Equations 

In this section we will outline the process to reduc e our 
correc t complete MHD equations to those found in IBelll 
(2004). As we stated at the end of section |2~2| the com- 
posite fluid formulation including the GCR inertia does 
not allow for MHD fluid instabilities to develop. How- 
ever, we can arrive at Bell's original system of equations 
if we do not properly account for the dynamics of the 
GCRs. 

We begin with the composite momentum equation 
(fT0|) , which recall included the inertial terms for the back- 
ground fluid, the excess electrons, and the GCRs. 



p— =— (V x B x B 
p dt 4ti- v ; 



(35) 



In addition to the composite momentum equation, we 
reintroduce the momentum equation for the GCRs. 

(7cr) n cr m p —^- = n cr e (~ u cr x B + E j (36) 
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Recall now equation (fTTj) . the electric field found above 
subject to the condition of negligible electron inertia, 
m e << mj. To simplify our subsequent calculations we 
will now consider the limit that n cr /nbi << 1. The mid- 
dle term in (|17|) is ~ n cr /n b i and can now be dropped. If 
we also neglect the Hall ter m (final term) from equation 
(fTTj) as done in [Bel (|2004f ), the expression for the elec- 
tric field simplifies to E = — (l/c)uu x B. Inserting this 
expression for E into equation (p5| the GCR momentum 
equation becomes the following. 

(7cr) n cr m p —j— = -j cr x B n cr eu bi x B (37) 

at c c 

Subtracting equation (|37[) from equation (|35[) we find the 
equation below. 



du du CT 

^~dt ^ cr ' ^cr^p 



(38) 



i(VxB)xB--j c 

47T C 



1 



x B H — n cr eu bi x B 



The LHS of this equation reduces to the inertial term for 
the background plasma and extra electron population, 
with the electron inertia neglected as before. 

du bl 1 
n bl m p —— =— (V x B) x B 
F cft 4-7T 

j cr x B + -n cr eu bl x B (39) 

c c 

This equation is identical to Bell's momentum equa- 
tion (equation ([2]) in this paper). Note that u&j in 
our formulation is the same as u in Bell's work. With 
E = — (l/ c)u&i x B, the induction equation also reduces 
to that of Bell (2004). 



dB 

~dt 



V x (u fcj x B) 



(40) 



Equations (|39|) and (|40)) represent two of the equations 
necessary for the development of the CRCD. We have 
found these equations from our composite formulation by 
removing the dynamics of the GCRs. Doing this has left 
out the GCR inertia from the above momentum equation 
but left information about the forces that act on them. 

In order to close the system of equations in this formal- 
ism, the dynamics of the GCRs must be a ccounted fo r 
in some way. This was the purpose of a in IBelll (2004). 
When a is set to 0, however, the information about the 
GCR inertia is lost. This would be equivalent to setting 
the LHS of equation (|3"Tf to zero, which implies that the 
final two terms on the RHS of equation (f39|) are also zero, 
removing the mechanism for driving the CRCD. 

We suggest that if the inertia of the GCRs is retained 
in the study of this problem it should balance the Lorentz 
forces on the background fluid responsible for driving the 
CRCD. In the event that the GCR inertia is neglected, 
the additional Lorentz forces are not consistent with the 
assumption that a — )• and should also be set to 0. 
This conclusion is more apparent given our method for 
considering the GCRs instead of the quantity a. 

This section does not imply that we should find Bell's 
instability in the calculations of Section l2.2.4l for the case 
where n cr /n b i << 1. The roots calculated in Section 
I2.2.4l were found for the case of the composite fluid where 



the GCR dynamics are included. In this section, we re- 
trieved Bell's equations by removing the GCR dynamics 
in addition to working in the limit of n cr /n b i << 1. 

3.3. Composite Formulation from Bell's Equations 

We can also retrieve our composite fluid momentum 
equation (fl"0|) from Bell's momentum equation (|2|). We 
begin with the momentum equation for the GCRs. 

du cr ( 1 \ 

(7cr) n cr m p —^- = n cr e I -u cr x B + E ) (41) 



dt 



As we have stated before, the electric field in Bell's for- 
mulation is equal to: 



E = -u x B 

c 



(42) 



Once again, u is the velocity of the background fluid. 
Inserting this expression for E into (|4ip the momentum 
equation for the GCRs becomes the following. 



du r 



1. 



x B 



1 



-rweu x B 



(43) 



dt c c 

Adding this to equation ^ we find that all of the Lorentz 
force terms cancel leaving only the momentum equation 
of ideal MHD for the composite fluid. This leaves no 
forces to drive an instability. By accounting for the in- 
ertia and dynamics of the GCRs we find that no MHD 
instability is excited by the presence of the GCRs. 

3.4. Additional Work by Previous Authors 

As we mentioned, work on this topic has been carried 
out by a number of authors in both the MHD and kinetic 
regimes. The purpose of this paper is to specifically ad- 
dress the previous linear MHD studie s of th e CRCD. In 
the MHD regime, iZirakashvili et all (|2008l ) have stud- 
ied the behavior of the instability. The authors confirm 
the existence of the instability proposed by Bell. This 
is not unexpected, however, because the MHD equations 
used to de scribe the problem are identical to those of 
[Belli ([200l . In this case, the electric current resulting 
from the drift of the GCRs exerts a magnetic Lorentz 
force on the background plasma which drives the insta- 
bility. In this paper, we propose a different method for 
analyzing this problem, within the confines of MHD. We 
have shown that if the inertia of the GCRs are treated 
on equal footing with that of the background plasma the 
relevant MHD equations take on a very different form. 
Most noticeable is the lack of a magnetic Lorentz force 
in the momentum equation which would drive the in- 
stability. We have shown that the total (magnetic and 
electric) Lorentz force responsible for driving the insta- 
bility is exactly balanced by the MHD inertial term of 
the GCRs. In the case of negligible GCR inertia, the 
total Lorentz force sums to zero leaving no net force to 
excite the instability. 

Kinetic and PIC studies have also been conducted with 
regards to the CRCD. These simulations, with and with- 
out feedback on the GCRs, show the possibility for mag- 
netic field amplification. A number of such simulations 
do not con sider parameters that include those stated in 
IBelll (120041) or those t hat re present realistic SNR condi- 
tions^ TGargat£eFal] (|2010l ) state that "The diverse pa- 
rameters chosen for the four sets of runs (A through D) 
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do not directly reflect typical physical conditions found 
in young SNR shock precursors for all the runs." Run 
A ment i oned here most closely resembles the conditions 
of I Belli (j2004l ) but does not consider a ny fee dback on 
the GCRs. The study of iNiemiec et al.1 (|2010t l consider 
a highly relativistic beam of ions, not a slowly drifting 
and nearly isotropic d istribution as we have don e here . 
The kinetic study by [Riquclme and Spitkovskvl (2009) 
also consider parameters very different from Bell. In 
that case, the GCRs stream at velocities > c/2 where as 
Bell considered the case of GCRs streaming at a speed 
of c/30. I n addit ion, the ratio of v a /u cr is ~ 10 ~ 4 in 
IBelll (12004ft while iRiauelme and Spitkovskvl (12009ft con- 
sider a ratio no less than 10~ 2 . O ur range of parameters 
includes those originally used by IBelll ([2004) as well as 
those relevant to SNRs. For this range we found no in- 
stability. Direct comparison between our work and those 
numerical simulations above is difficult as the conditions 
in those papers consider GCR flow speeds approaching 
c. The use of relativistic flow speeds is not valid in our 
non-relativistic MHD approximation. This is why we 
have chosen to compare our results directly with those 
of the original MHD situation considered by Bell. 

4. DISCUSSION AND CONCLUSIONS 

In this paper we have considered the problem of MHD 
fluid stability in a plasma permeated by a fluid of stream- 
ing GCRs. This work was carried out in the framework 
of MHD which is valid in typical astrophysical situations. 
When considering the effects of the GCRs on the back- 
ground plasma it is important to consider how the GCRs 
are treated in the problem. Here we have discussed two 
different methods of approaching the problem which lead 
to very different results. 

If the electric current created by the streaming GCRs 
is simply added to the dynamical equations for the back- 
ground plasma without considering their own inertia, as 
it appears to have previously been done, there is a pos- 
sibility for an instability to develop if condition (1321) is 
satisfied. The instability greatly enhances the previously 
zero transverse component of the magnetic field. Physi- 
cally, the instability is the result of the perturbed mag- 
netic field being stretched out by the Lorentz force be- 
tween the spiraling field and the neutralizing return cur- 
rent in the background plasma. The — j cr X B term in ([2]) 
acts as a driving force to what would otherwise be simple 
Alfven waves. For a brief discussion of the physics of the 
instability as well as an illustrative figure of the relevant 
fo rces see Fig . 1 in[Zirakash vih' et al.l (|2008ft . In the case 
of IBelll (|2004ft , j cr is constant and free from any feedback 
allowing it to continually drive the return current and 
the instability. 

Here, we have chosen to re-examine the original MHD 
analysis while explicitly including the inertia and dynam- 
ics of the GCRs which has been neglected in previous 
MHD studies. In this paper, we have shown that the 
method for including the GCRs has a dir ect bearing on 
the cr eation of an instability. As stated in INiemiec et al.l 
(2008) , "Previously published MHD simulations have as- 
sumed a constant cosmic-ray current and no energy or 
momentum flux in the cosmic rays, which excludes a 
back-reaction of the generated magnetic field on cosmic 
rays, and thus the saturation of the field amplitude is ar- 
tificially suppressed." These authors found no evidence 



of the CRCD instability and only a moderate amplifica- 
tion of the magnetic field (SB/B ~ 1) due to turbulent 
motions in the plasma. We have shown in this paper 
that failure to consider the dynamics of the GCRs leads 
to excess forces on the system (the j cr x B term in ([2])) 
which stretch out the magnetic field lines. This force 
is balanced by the GCR inertial in our composite fluid 
formulation. 

This paper highlights the importance of how the GCRs 
are accounted for in this problem. By considering the dy- 
namics of the GCRs (and thus their inertia) along with 
that of the background fluid we find that there is no 
such MHD fluid instability. By working with the com- 
posite fluid, the MHD momentum equation contains the 
total current, j, in the system rather than only the re- 
turn current — j cr . In the case of an equilibrium uniform 
magnetic field j = so there is no driving term as in 
Bell's momentum equation ((2|) to excite an instability. 
Physically, the lack of any net current means that there 
is no net Lorentz force to stretch out the perturbed mag- 
netic field. We have also shown that the CRCD is the 
result of neglecting the GCR inertia in the analysis of 
this problem. 

The lack of an instability in our analysis calls into ques- 
tion our neglect of the kinetic resonant streaming insta- 
bility due to its slow growth rate compared to that of 
the non-resonant instability. This concern will be ad- 
dressed through a numerical study of this problem in a 
subsequent paper. 

The conclusion of t his paper implies that th e observed 
large magnetic field ([Berezhko fc Volkl 12003ft in SNRs 
should be found in some process other than cosmic-ray- 
driven amplification upstream of the blast wave. One 
such process is the downstream a mplification caused by 
pre-existing upstream turbul ence (jGiacalone and Jokipiil 
2007; Bcr esnvak etHl 12009ft . In this scenario, the ac- 
celeration to energies ~ 3 x 10 15 eV could occur at the 
quasi-perpendicular parts of the blast wave0 

This work was supported, in part, by NASA under 
grants NNX08AH55G and NNX10AF24G. JRJ is grate- 
ful to A. R. Bell and M. A. Lee for helpful discussions 
which improved this paper. 
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APPENDIX 

A. DERIVATION OF COMPOSITE FLUID MOMENTUM EQUATION 

The ith component of the four momentum equations are shown below in index notation. In order the equations are 
for the background ions, background electrons, GCRs and the neutralizing electrons. 











— + m v n hl u hi ^ + ^— [m p n b i (u bi ,iUj + uu,jUi - UiUj)} = n bi e 



dt 



d d , 

T^m p nbiUbe,i + -q^- [m e nu [Ub e ,iUj + Ub e ,jU l — u%Uj)\ — —n b ie 







Ei H — eijkUbe.jBk 



d ., 



dt 



dt 



TflpTicr^ex ,i 



+ 



8x4 



\jTlpTl cr {u^Xfi^j He 



UiUj)} 



Ei -\~ €ijk'Uex,j Ek 



Summing equations (|A1[) - (|A4[) we arrive at the following equation for the composite fluid. 

^/i ^ . . 1 . _ 
^{pUij + -^-(puiUj) = -ei jk JjB k 



dt 



(Al) 
(A2) 
(A3) 
(A4) 

(A5) 



Recall the definitions of p and B from section 12.21 We then expand the derivatives on the LHS of the above equation. 



dm Ui 
p -dt +fMj d^ 



dp dpuj 

Ui— +Ui 



1 



at 



Factoring the third and fourth terms we see that these terms cancel by mass conservation. 



dui 



8x4 



dp_ dpu-j 
dt dxj 



(A6) 



(A7) 



We now substitute for the total current j by Ampere's law to arrive at our momentum equation for the composite 
fluid, equation (flC)]) above. 



^ = I( V xB)xB 



dt 4?r 
B. CALCULATION OF {y CR ) 

The factor j cr is an average Lorentz factor for the entire GCR population and can be calculated as shown below 

III lcrf(p)dp 



(A8) 



(lev) 



(Bl) 



If we assume the GCR distribution function f(p) to be a power law, ~ p 1 , we find the following form for (7 cr ). 



(7cr> 



pp-r 17-E-V 



-r 

-,1/2 



+ 1 



dp 



(B2) 
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The limits of integration p\ and pi are the minimum and maximum momentum of the CRs in the distribution function 
f(p). For T = 4 as in the case of ideal DSA at a strong shock this integral can be solved analytically. 



(lev) 







3/2 










\m p cj 





(B3) 



